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Abstract  The  paper  surveys  some  of  the  results  related  to  approach  of 
hierarchic  modeling  of  the  plate  problems.  The  main  ideas  are  explained  and 
illustrated  by  numerical  examples. 
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1.  Introduction 


The  problem  of  partial  differential  equations  on  a  special  domain  w, 
typically  the  "thin"  domain  is  usually  simplified  by  various  dimensional 
reduction  techniques.  The  aim  is  to  approximately  solve  the  3  dimensional 
problem  by  a  two  dimensional  formulation.  This  approach  is  widely  applied  in 
the  connections  with  plates  and  shells.  Today  many  plates  and  shells  theories 
(models)  are  used  in  practice. 

The  principles  of  derivation  of  these  models  can  be  divided  into  3 
groups . 

a)  Physical  derivation.  Here  various  "a  priori"  assumptions  of  geo¬ 
metrical  nature  are  made  together  with  some  additional  assumptions  involving 
Importance  of  certain  stress  components.  A  typical  example  is  the  Klrchhoff 
hypothesis  Ill.  This  approach  was  generalized  by  combination  with  various 
heuristic  approaches  to  treat  the  anisotropic  and  laminated  plates  and  shells. 
We  mention  [2],  13],  [4]  as  an  example  of  this  type  of  derivation  (see  also 
references  there). 

b)  Asymptotic  analysis  Here  the  main  idea  is  to  use  the  power  series 
expansion  (in  the  thickness)  of  the  solution.  It  was  analyzed  with  various 
levels  of  rigor.  We  mention  here  [5-12]  and  references  mentioned  there  as 
typical  for  this  approach. 

Many  results  were  devoted  to  the  studies  of  the  convergence  (when  the 
thickness  d — >0)  and  the  accuracy  of  the  solution  and  the  arguments  why 
these  models  are  relevant.  We  mention  here  for  example  [4],  [11-16]  and 
reference. s  there. 

The  appioach  of  this  type  lead  to  the  formulations  of  well  known  plates 
models,  as  for  example,  the  Relssner-Mindlin  model  initiated  in  [8]  [17]  and 
which  is  a  special  case  of  the  Naghdi  plate  theories  [11]. 
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Various  plate  models  require  various  types  of  boundary  conditions;  more 
exactly  some  models  allow  to  distinguish  between  certain  boundary  conditions, 
while  others  do  not.  As  an  example  we  mention  that  the  Reissner-Mindlin  model 
allows  to  distinguish  between  a  soft  and  hard  simple  support,  while  the 
Kirchhoff  model  does  not.  These  various  boundary  conditions  could  lead  to 
very  different  properties  of  the  solution.  As  an  example  we  mention  the 
paradox  for  the  hard  simple  support  [18-22]  which  says  that  the  solution  of 

4) 

the  hard  sample  support  on  a  regular  n-gone  inscribed  in  the  unit  circle  does 
not  converge  to  the  solutions  on  the  circular  domain  as  n — >».  In 
contrast,  for  the  soft  simple  support,  this  paradox  is  not  present.  The  hard 
support  paradox  also  occurs  in  the  3  dimensional  theory  [21]. 

The  properties  of  solution  of  the  Reissner-Mindlin  model  especially  Its 
boundary  layer  behavior  for  smooth  domain  was  studied  rigorously  in  [23] , 

[241.  In  [25]  the  boundary  layer  for  unsmooth  domain  was  heurlstically 
analyzed.  The  behavior  of  the  solution  in  the  neighborhood  of  the  corners  of 
the  domain  was  analyzed  in  [26]  [27].  For  the  survey  of  various  theories  for 
laminated  plates  we  refer  also  to  [4]. 

c)  Numerical,  hierarchical  approaches.  This  approach  is  a  modern  one, 
and  typical  for  the  computer  oriented  procedures.  It  is  stemming  from 
general  principles  of  numerical  methods,  adaptive  approaches  and  a-posteriori 
error  estimations.  Here  no  particular  model  is  a-priori  preferred.  This 
approach  constructs  adaptively  an  optimal  numerical  method  for  solving  the 
three  dimensional  problem  which  leads  to  the  accurate  computation  of  data  of 
interest. 

In  the  next  s^cHons  we  will  elaborate  on  this  approach.  For  the  sake 
of  concreteness  we  will  address  these  questions  on  the  specific  model  problems 
although  the  results  hold  in  a  general  setting. 
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2)  The  three  dimensional  formulation  of  the  plate  problem 


Let  «eR  be  a  bounded  domain  with  the  boundary  T  and  £1  =  {x  = 
3  di  d 

(xltx2,x3)  €  R  |  (xltx2)  €  u,  -|<x3*|).  Further  we  let 


, x2) €  r,  ~  2  ^ X3 < 


S  =  jxe  R3|  (Xj, 

R±  =  jxe  R3|  (Xj.x^  €  «,  x3  =  ±dj. 


We  will  refer  to  the  solution  of  the  three  dimensional  elasticity  problem  for 
an  isotropic  homogeneous  material  when  the  equal  normal  load  |q(x)  is  acting 
on  R+,  as  the  exact  solution  of  the  plate  problem.  The  solution  u  = 
(Uj.iu.u^  is  the  minimizer  of  the  total  energy 

(2.1)  GA(u)  =  ^(u)  -  Q(u) 


1  3 

over  the  set  of  functions  X(ft)c  (H  (Q))  which  satisfy  certain  constraints 

on  S.  The  boundary  conditions  of  the  plate  problems  are  uniquely 

characterized  by  this  set  X(ft). 

A 

In  (2.1)  we  denoted  by  &  (u),  the  strain  energy  based  on  the  Hooke’s 
compliance  matrix  A  xor  the  plate  material. 

The  problem  has  been  formulated  for  simplicity  only  for  the  homogeneous 
boundary  conditions.  For  concreteness  and  simplicity  we  will  deal  with  the 
square  plate  which  we  denote  by  ft  resp.  u  =  (x1tx_|  |x, I  <0.5,  i  =  1,2} 

S  S  1  Z  1 

and 

Tj  *  <X!.X2I  xj  =  0-5,  I x2 1  <0.5) 
r2  =  (Xj , x2 1  Xj  <  0.5,  x2  =  0.5) 
r3  =  (Xj , x2 1  x  1  =  -0.5,  i x2 1  <0.5) 
r4  =  (Xj.x^  I Xj |  <  0.5,  x2  =  -0.5) 

S1  =  rix  (-d/2,  d/2),  i  =  1,2,  •••,4. 
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The  constrain  X(Q)  defines  various  boundary  conditions.  Let  us  show 


some  typical  ones  which  could  be  imposed  for  example  on  x  (-d/2,  d/2) 

a)  clamped 

b)  free 

c)  hard  simple  support: 

d)  soft  simple  support: 


U1  =  U2  =  U3  =  0  °n  S1 


no  constrain  on  S, 


u2  =  u3  =  °  on  Sj 


u^  =  0  on 


d/2 


e)  very  soft  simple  support: 


u^fXj.x^.x^j)  dx^  =  0,  (Xj.x^sfj 


-d/2 


f)  illegal  simple  support:  u^Cx^x^O)  =  0,  (x^x^efj  . 


The  last  type  of  the  simple  support  is  called  " illegal" ,  because  its 

solution  is  the  same  as  for  the  free  boundary  condition.  Mathematically, 

1  3 

this  follows  from  the  fact  that  functions  of  (H  (£1))  do  not  have  trace  on 
the  one  dimensional  line  =  {x^.x^.x^  I  xj>x2€  x3  = 

Physically  it  follows  from  the  observation  that  the  displacement  under 
load  concentrated  on  a  line  is  infinite  (Bousinesque  solution). 

The  solution  of  the  three  dimensional  problem  with  finite  strain  energy 
exists  and  is  uniquely  determined  (except  of  the  case  of  free  boundary 
conditions  on  T  when  the  usual  additional  conditions  have  to  be  imposed). 

The  aims  of  the  computation  can  be  different.  For  example  the  interest 
can  be  in  the  bending  moments  and  shear  forces  or  average  displacements,  etc 
Further  the  interest  could  be  only  in  these  values  inside  the  domain  w  or 
at  the  boundary  T  or  in  the  neighborhood  of  T.  We  can  be  interested  in 
the  stresses  or  in  the  stress  Intensity  functions  along  the  edges  and  stress 
intensity  factors  in  the  vertices  etc. 
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It  is  obvious  that  any  imposed  boundary  condition  is  an  idealization. 
Hence  one  cam  study  how  the  results  are  influenced  bv  the  uncertainty  of  the 
idealization  of  the  boundary  conditions.  We  will  not  address  this  problem  in 
details  here.  Instead  we  show  am  example  of  the  modeling  of  the  clamped  in 
boundary  condition  for  the  plate  in  the  case  when  the  solution  is  independent 
of  the  variable  x^. 

Let  us  consider  the  plate  of  thickness  d  =  1  and  span  L  =  20  with 
the  modeled  clamped  end  as  shown  in  Fig.  2.1.  We  will  assume  that  the 

7 

material  is  homogeneous,  its  modulus  of  elasticity  E  =  3. 10  and  Poisson 
ratio  v  =  0.3.  Boundary  condition  e)  is  the  spring  condition  with  the  spring 

g 

constamt  C  =  10  .  Case  f  in  the  Fig.  2.1  is  the  strength  of  material 
model  based  on  the  K.irci*hoff  theory.  In  this  case  the  bending  moments  in  the 
middle  of  the  plate  is  M  *  ^  202q  =  q.  Normalizing  the  results  by  the 

moment  of  the  model  d)  we  give  the  error  in  the  bending  moment  M  in  the 
middle  of  the  plate  in  the  Table  2.1. 


Figure  2.1  The  scheme  of  various  boundary  conditions 
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Table  2. 1  The  error  of  various  models  of  boundary  conditions 


Case 

Error  in  M 

a 

+  8.2*/. 

b 

+  9.5% 

c 

+  8.0% 

d 

0% 

e 

+19.9% 

f 

-  0.2% 

Table  2.1  clearly  shows  the  influence  of  the  uncertainties  in  the  boundary 
condition  on  the  bendings  moment.  The  influence  on  the  stresses  in  the 
neighborhood  of  the  boundary  condition  is  much  larger.  Tor  more  see  [28]. 


3.  The  principles  of  the  numerical  hierarchical  modeling 
Let  us  assume  the  following  form  of  the  displacement 

ni 


(3.1) 


ui(xrx2,x3)  =  i>i,j<xrV  ^i,  j(^d  ] ’ 
j=0 


i  =  1,2,3 


where  the  functions  <p.  .  (tj),  -1<tj<1,  may  in  general  depend  on  d  and 

*  *  J 

n  =  (n^  n2,n^),  n^O,  1  =  1,2,3  are  integers.  Because  of  our  assumptions 

about  the  symmetry  of  the  loading  we  can  assume  that  <p.  .(tj)  are  anti 

*  *  J 

-symmetric  functions  for  i  =  1,2  and  symmetric  for  j  =  3.  Then  u,  =  0  for 

*  »  J 

1  =  1,2  and  j  even  and  u_  =  0  for  J  odd.  We  will  assume  that 

3.  J 

1  3 

u €  (H  (Q))  nK(fl)  and  n  can  be  different  in  different  regions  of  w. 

The  following  questions  arise 

a)  how  to  select  the  functions  <p.  (tj) 

1  »  J 

b)  how  to  select  n;  uniform  or  with  different  values  in  different 
regions  of  w 


6 


c)  how  to  derive  the  differential  equations  or  variational  formulations 
to  determine  the  functions  u,  . . 

d)  how  to  assess  the  accuracy  of  the  solution  with  respect  to  the  data 
of  interest  (when  compared  with  the  data  of  exact  3  dimensional 
solution),  how  to  design  an  adaptive  procedure. 

Specifying  the  functions  <p.  (tj )  and  the  principles  of  derivation  of 

1  #  J 

the  differential  equations  for  a  sequence  n^  =  (n|^,n|^,nj^),  1  =  1,2, 
and  any  thickness  d  we  construct  a  family  of  n-models.  This  family  will  be 
called  hierarchical  if  it  satisfies  the  following  conditions. 

i)  Let  ud  and  und  be  the  solution  of  the  3  dimensional  problem  of  the 
plate  of  thickness  d  and  the  solution  of  the  n-model  respectively.  Then 

(3.2)  l|ud  “  “nd^^d® — *°  as  d — *° 
holds. 

Mostly  but  not  exclusively  the  energy  norm  ||  •  ||  in  (3.2)  is  considered. 

il)  if  the  solution  u^  is  sufficiently  smooth  uniformly  with  respect  to 
d  then 

(3.3)  ,ud“un  dl/|udl  SC(n1)d*(n,) 

where  a(ni+1)  >  atnj) 

iii)  for  any  fixed  d  we  have 

(3.4)  «ud-un  dll/«udll — >0  as  i — 

and  preferably  but  not  necessarily  ||u. -u  ,||£[|u-u  ,|| . 

d  riyd  ni+id 

We  remark  that  the  notion  of  hierarchy  depends  strongly  on  the  norm  ||  •  ||  . 
We  can  also  define  the  hierarchic  family  by  replacing  the  norm  ||  •  ||  by  another 
definition  of  accuracy,  which  expresses  the  goal  of  the  analysis. 
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T 


4.  The  basic  hierarchy 

Let  us  assume  that  functions  <p,  .(■»))  in  (3.1)  are  independent  of  d. 

*  »  J 

Then  in  the  spirit  of  the  requirement  (3.3).  we  can  ask  which  <f>.  ’s  lead 

1  *  J 

to  that  highest  a.  The  answer  is:  <p.  .  have  to  be  polynomial  of  degree  j. 

*  »  J 

For  the  arguments  leading  to  this  conclusion  we  refer  to  [29],  [30],  [31]. 
Hence  (3.1)  can  be  written  as 


(4. 1) 


ni 


ui(x1,x2,x3) 


j=0 


X2} 


1,2,3 


and  as  before  u  =0  for  i  *  1,2  and  j  even  and  u„  =  0  for  j 
1  ^  »  J 

odd.  The  main  models  we  can  therefore  consider  are:  n  =  (1,1,0),  (1,1,2), 
(3,3,2),  (3,3,4)  etc. 

Let  us  now  address  the  question  of  deriving  the  differential  equation 

for  u  We  define  u,  .  as  the  minimizers  of  the  potential  energy 

*  >  J  *  *  J 


(4.2) 


GB(u)  =  gB(u)  -Q(u) 


over  all  u  of  the  form  (4.1)  belonging  to  X(n)c  (H^Q))3  (see  also  (2.1)) 
and  B  boing  a  Hooke’s  compliance  matrix  which  has  not  to  be  identical  with 
the  matrix  A  of  the  plate  material.  We  have 


(4.2) 


(l-i>)E  i>E  vE 

(l+i>)(l-2w)  (l+v>)d-2i>)  (1+v)  (l-2v) 

( 1  -v)E  vE 

(l+v)(l-2i>)  (l+i>Hl-2i>) 


A  = 


Symmetric 


( l-i>)E 
(l+w)(l-2v) 


0 

0 

0 

E 

1+v 


0 

0 

0 

0 

E 

1+v 


0 

0 

0 

0 

0 

E 

l+i> 
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As  usual,  we  denote  by  E,  the  modulus  of  elasticity,  and  by  v,  the 
Poisson  ratio. 


We  have 

Theorem  4. 1  The  model  (1,1,0)  with 


(4.4) 


VI  s  u  2X3 

1  1,1  d 

2x3 

u  =  u  — - 

2  2,1  d 


U3=U3,0  and 
B  =  A 

is  for  v*0  not  the  member  of  the  hierarchic  family  (it  violates  condition 
(3.2)). 

For  more  see  [321,  [33]. 

Let 

E  vE  _  _ 


(4.5) 


iV  1-v2 


1C 

-V 


0 

E 


1-v 


Symmetric 


E 

1+v 


Ek 

l+v 


Ek 

l+v 


where  0 < k < 1  arbitrary.  Then  we  have 


Theorem  4.2  The  model  (1,1,0)  is  the  member  of  the  hierarchic  family  when 
the  matrix  B  *  R  defined  by  (4.5)  is  used.  Then  for  smooth  solution 
II ud  "  UndM ud^  *  ^  <**1.  Model  (1,1,0)  with  matrix  B  =  R  is  the 

well  known  Reissner-Mindlln  model  (see  [32],  [33]). 
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The  coefficient  k  in  (4.5)  is  the  shear  factor  with  different  values 


recommended  in  the  literature.  The  typical  recommendation  is  k  =  5/6  or 

0.833 < k < 0.870  as  recommended  in  [34]. 

Let  us  consider  the  example  of  square  plate  Q  (see  Section  2.1) 

s 

7 

uniformly  loaded  with  simple  soft  support,  and  d  *  0.01.  Assume  E  -  10 
and  v  *  0.3.  Further  by  ?{k)  we  denote 


(4.6) 


SU) 


1/2 


|g(»c)  -  8(3  dim) | 
|g(3  dim)  | 


where  by  8( 3  dim)  we  denoted  the  strain  energy  of  the  3  dimensional  solution 
and  by  £(*)  the  strain  energy  of  the  (1,1,0)  model  with  B  =  R  given  by 
(4.5). 

In  the  Table  4.1  we  show  the  energy  8(k),  8(3  dim)  for  various  k  and 
the  error  £(*). 

Table  4.1  Energy  S(k)  as  function  of  ic  and  the  error  £. 


K 

g(K)  x  103 

C/. 

1 

0 . 234563 

2.69 

0.91 

0 . 234674 

2.26 

0.87 

0 . 234729 

0.41 

0.8453 

0 . 234765 

1.16 

5/6 

0 . 234783 

1.46 

KIRCHHOFF 

0 . 232392 

9.98 

3  DIM 

0 . 234733 

— 

In  the  Table  4.1  we  also  reported  the  solution  of  Kirchhoff  model.  We  see 
that  k  *  0.87  leads  here  to  the  best  result  and  the  Kirchhoff  model  has  a 
large  error  in  the  energy  norm. 

Remark  The  solution  of  the  3  dimensional  formulation  is  not  known  In  the 
analytic  form.  This  solution  was  numerically  computed  by  the  hp-version  of 
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FEM  with  refined  meshes  and  error  estimation  which  guarantee  that  the  data 
reported  here  are  adequately  accurate,  i.e.,  the  error  does  not  influence  our 
conclusions. 

As  an  illustration  of  Theorem  4.2  consider  the  problem  of  clamped-in 
plate,  uniformly  loaded.  In  the  table  4.2  we  show  £(0.87)  as  function  of  d 

Table  4.2.  The  energy  6(k)  and  the  error  £  as  function  of  d. 


d 

3  DIM 

RM,  *c=0.S7 

B 

0.242115(-6) 

0.245521 (-6) 

11.8*/. 

0. 1491 15 ( —4 ) 

0. 149256 (-4) 

3.07% 

69 

0.232489 (-3) 

0.232524 (-3) 

1 . 2% 

We  clearly  see  the  convergence  rate  a  ■  1. 

E 

The  shear  factor  influences  only  the  two  terms  in  the  lower  right 

corner  of  the  matrix  A.  Hence  we  could  consider  also  for  higher  models  the 

* 

matrix  A  with  the  shear  factor  influencing  these  two  values. 

We  have 

Theorem  4.3  The  models  (n,n,n+l),  n  =  1,2,***,  using  the  matrix  B  =  R  do 
not  belong  to  the  hierarchic  family  because  condition  (3.2)  does  not  hold. 

Theorem  4.4  The  sequence  of  the  models  (n,n,n+l),  n  =  1,2,* ••  is  hierarchic 
for  B  =  A  and  is  not  hierarchic  when  B  ■  A  . 

Table  4.3  which  addresses  the  problem  of  the  square  plate  with  soft  support 


illustrates  Theorem  4.4. 


Table  4.3  Energy  8  for  different  models  and  the  error  £ 


K 

Model 

8  •  103 

C/. 

B 

1 

(1,1,2) 

0 . 234528 

2.95 

A 

1 

(3,3.4) 

0 . 234731 

0.29 

A 

1 

(5,5,6) 

0.234732 

0.01 

A 

* 

0.87 

(1,1,2) 

0 . 234693 

1.3 

A 

0.87 

(3,3,4) 

0.234913 

2.76 

A 

* 

0.87 

(5,5,6) 

0.234914 

2.78 

A 

3  dim 

0 . 234733 

Remark  In  all  examples  we  mention  the  exact  solutions  of  the  differential 
equations  of  the  models.  They  were  solved  numerically  with  the  error  control 
so  that  the  reported  data  are  exact  in  the  range  we  used  for  our  conclusions. 


5)  Capturing  the  boundary  layer 

As  has  been  said  above,  we  compare  the  solution  of  the  n-model  with  the 

solution  of  the  three  dimensional  problem.  The  accuracy  then  has  to  be 

assessed  relatively  to  the  data  of  interest. 

Let  us  be  interested  in  the  behavior  and  accuracy  of  the  bending  moments 

and  shear  forces  in  the  neighborhood  of  the  boundary. 

Consider  the  problem  of  unit  square  plate  ft  (d  =  0.01)  and  compare  the 

s 

values  of  the  twist  moment  computed  from  the  3  dimensional 

formulation,  RM  model  (i.e.  (1,1,0)  model  with  B  =  R  and  k  =  0.87)  and 
(1,1,2)  model  with  B  =  A,  k  *  1. 


12 


Table  5.1a  Twist  moment  M12(0.4,x2)  for  the  hard  and  soft  simple  support. 


0.005 


0.030 


Table  5.1b  Twist  moment  M12(0.5,x2)  for  hard  support. 


3  DIM 


0.002 

0.010 

0.017 

0.031 

0.032 

0.032 

0.032 

0.033 

0.033 


0.033 

0.035 
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We  see  clearly  that  the  boundary  layer  is  present  for  the  soft  support 


and  Is  practically  not  existent  for  the  hard  support.  Let  us  remark  that 
M12(0-5,X2)  =  0  for  the  soft  support. 

In  the  tables  5.2a,  b  we  show  the  analogous  results  for  the  shear  forces 
Q3I  and  Q32- 


Table  5.2a.  The  shear  forces  Q31(0.4,X2)  for  the  hard  and  soft  support 


HARD 

SOFT 

X2 

3  DIM 

RM 

(1,1,2) 

3  DIM 

RM 

(1,1,2) 

0 

0.246 

0.246 

0.246 

0.245 

0.246 

0.246 

0.02368 

0.245 

0.245 

0.245 

0.245 

0.245 

0.245 

0.11842 

0.235 

0.234 

0.234 

0.235 

0.234 

0.234 

0.21316 

0.208 

0.208 

0.208 

0.208 

0.208 

0.208 

0.45000 

0.052 

0.053 

0.053 

0.053 

0.056 

0.055 

0.48079 

0.020 

0.021 

0.021 

-0.005 

0.001 

0.006 

0.49026 

0.010 

0.011 

0.011 

-0.419 

-0.403 

-0.343 

0.49500 

0.005 

0.005 

0.005 

-1.921 

-1.912 

-1.822 

0.49713 

0.003 

0.003 

0.003 

-3.775 

-3.817 

-3.824 

0.49903 

0.001 

0.001 

0.001 

-6.974 

-7.047 

-7.379 

0.49950 

0.001 

0.000 

0.000 

-8.167 

-8.214 

-8,695 

0.50000 

0 

0.000 

0.000 

-9.718 

-9.639 

-10.323 
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Table  5.2b  The  shear  force  Q_.  (0.5,x_)  and  Q„,(0.5,x_)  for  the  soft 

Jx  £.  c. 

support 


Q33 

l(0.5.x2) 

Q32(0.5,x, 

>) 

X2 

3  DIM 

RM 

(1.1,2) 

3  DIM 

RM 

(1,1,2) 

0 

0.420 

0.420 

0.420 

0 

0 

0 

0.02368 

0.420 

0.420 

0.420 

-0.643 

-63801 

-0.685 

0.11842 

0.408 

0.408 

0.408 

-3.199 

-3. 173 

-3.399 

0.21316 

0.374 

0.374 

0.374 

-5.663 

-5.617 

-6.018 

0.45000 

0. 133 

0.133 

0. 135 

-10.451 

-10.363 

-11.005 

0.48079 

-0.133 

-0. 125 

-0.105 

-10.635 

-10.559 

-11.332 

0.49026 

-0.813 

-0.807 

-0.756 

-9.905 

-9.844 

-10.651 

0.49500 

-2.175 

-2.188 

-2.157 

-8.600 

-8.505 

-9.251 

0.49713 

-3.514 

-3.573 

-3.645 

-7.706 

-7.569 

-8.221 

0.49903 

-5.462 

-5.499 

-5.815 

-6.954 

-6.881 

-7.368 

0.49950 

-6.069 

-6.073 

-6.485 

-6.838 

-6.727 

-7.252 

0.50000 

-6.793 

-6.678 

-7.221 

-6.793 

-6.678 

-7.221 

Let  us  analyse  now  the  strength  of  the  boundary  layer.  To  this  end  let 


us  define  the  functions  031(x1>x2). 


exp 


^-031(x1(x2)  (0.5-x2)/dj 


°31(xrx2*  ~  ^31^X1,X2^ 
^(Xj.0.5) 


Because  the  hard  support  solution  has  no  boundary  layer,  we  use  it  as  a  "base" 
smooth  function.  Functions  03j  characterizes  the  strength  of  the  boundary 
layer.  Table  5.3  shows  the  values  of  03^(x^,x2)  for  x^  =  0.4,  0.5. 
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Table  5.3  The  strength  function  of  the  boundary  layer  for  the 

soft  support 


Xj  =  0.4 

x^  =  0.5 

X2 

3  DIM 

RM 

(1,1,2) 

3  DIM 

RM 

(1,1,2) 

0.49263 

3.22 

3.25 

3.43 

2.17 

2.21 

2  37 

0.49500 

3.24 

3.24 

3.47 

2.27 

2.23 

2.41 

0.49666 

3.28 

3.24 

3.47 

2.28 

2.20 

2.39 

0.49878 

3.40 

3.23 

3.47 

2.23 

2.04 

2.24 

0.49950 

3.51 

3.23 

3.47 

2.22 

1.91 

2. 12 

0.49989 

3.62 

3.24 

3.47 

2.22 

1.83 

2.12 

0.49991 

3.63 

3.24 

3.47 

2.21 

1.83 

2.03 

0.49994 

3.63 

3.23 

3.47 

2.21 

1.83 

2.03 

It  has  been  shown  In  [23]  [24]  that  for  the  smooth  boundary  P31(xlfx2) 

=  V12k  =  3.23  (in  our  case)  which  can  be  used  if  Xj  *  0.5.  In  [25]  is 
suggested  that  £^(0.5,  x ~  V6ti  =  2.28.  We  see  that  the  value  3.23  very 
well  describes  the  character  of  the  boundary  layer  while  the  exact  theoretical 
strength  for  x^  =0.5  is  yet  an  open  question.  We  also  see  that  the 
boundary  layer  of  the  3  dimensional  solution  as  well  as  of  the  model  (1,1,2) 
is  stronger  than  of  the  RM  model.  We  will  return  to  this  question  later  in 
this  section. 

As  d — >0  then  the  difference  between  hard  and  smooth  support  dis¬ 
appears  when  measured  in  the  energy  norm.  It  is  not  the  case  for  the  data  at 
the  boundary.  We  can  expect  for  example  that  dQ3^(0.5,  0.5)  converges  as 
d — >0  to  some  value  (although  not  yet  proven). 

In  the  table  5.4  we  show  the  values  of  Q31(0.5,  0.5)  and  dQ31 (0.5,0. 5) 
computed  from  RM  model  for  k  =  5/6. 
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Table  5.4  The  values  of  <2^(0. 5,  0.5)  for  RM  model  as  function  of  d  for 
the  soft  simple  support. 


d 

Q31(0.5,  0.5) 

dQ31(0.5,  0.5) 

0.025 

-2.68 

-0.0671 

0.01 

-6.54 

-0.0654 

We  can  expect  that  for  fixed  d,  Q^(0.5,  0.5)  is  proportional  to  Vk. 
In  Table  5.5  we  show  some  values  for  d  =  0.01. 


Table  5.5  Q^CO.S,  0.5)  as  function  of  k  for  the  RM  model 


K 

Q31(0.5,  0.5) 

k~U2  Q31(0.5,  0.5) 

1 

-7.15 

-7.15 

0.91 

-6.83 

-7.16 

0.87 

-6.68 

-7.16 

5/6 

-6.54 

-7.16 

So  far  we  have  seen  that  the  RM  model  describes  reasonably  well  the 
moments  although  the  strength  of  the  boundary  layer  is  more  than  10%  off. 

Our  discussion  was  for  the  soft  simple  support.  For  other  boundary 
conditions  the  situation  can  be  different.  As  an  example  we  consider  the 
problem  of  the  square  plate  of  thickness  d  =  0.01  with  all  four  sides 
clamped-in.  In  the  Table  5.6  we  report  the  values  of  the  moments  M^tx^.O) 
and  M22(xlt0). 
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Table  5.6  The  moment  M^fx^O)  and  M22(Xj,0)  for  various  models 


MjjtXj.O) 

M22l*r0) 

X1 

3  DIM 

RM 

(1,1,2) 

3  DIM 

RM 

(1,1,2) 

0 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

-0.0229 

0.0000 

-0.0157 

-0.0157 

-0.0157 

-0.0157 

-0.0163 

-0.0163 

0.4000 

+0.0164 

+0.0163 

+0.0164 

+0.0026 

+0.0027 

+0.0026 

0.4900 

0.0470 

0.0470 

0.0470 

0.0141 

0.0141 

0.0141 

0.4930 

0.0483 

0.0483 

0.0483 

0.0145 

0.0144 

0.0144 

0.4990 

0.0509 

0.0509 

0.0509 

0.0168 

0.0152 

0.0170 

0.4993 

0.0510 

0.0510 

0.0510 

0.0176 

0.0153 

0.0179 

0.4999 

0.0512 

0.0512 

0.0512 

0.0207 

0.0153 

0,0212 

0.5000 

0.0513 

0.0513 

0.0513 

0.0220 

0.0153 

0.0220 

We  see  that  the  error  in  M^Cx^O)  is  negligible  for  RM  and  the  (1,1,2) 
model  and  no  boundary  layer  is  present.  The  moment  M22^X1 , °)  computed  from 
the  RM  model  has  no  boundary  layer  either,  but  has  error  30%.  In  contrast, 

3  dimensional  solution  shows  very  strong  boundary  layer. 

We  will  analyze  the  strength  of  the  boundary  layer  as  before.  Let  us 
define  822(xltx2) 

.xp[-g22(x1.x2)[^i^]]  . 

I  *^2^X1*X2^  “ M22^x1’x2^  I 

ll^(0.S.x2)-M22(0.5.x2)1  ' 

Here  we  did  use  •S'  Xj,x2)  as  the  smooth  extension  of  M??. 
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Table  5.7  The  strength  function  022^X1’^  °*  ^ovindary  layer  for  the 

clamped  plate. 


X1 

3  DIM 

(1,1,2) 

0.4990 

15.43 

13.86 

0.4996 

18.17 

13.41 

0.4999 

23.85 

13.32 

0.49993 

24.13 

13.32 

0.49999 

24.78 

13.30 

We  see  that  the  3  dimensional  solution  shows  very  strong  boundary  layer 
for  the  moment  M^.  The  model  (1.1.2)  shows  a  relatively  strong  boundary 
layer  and  reasonable  accuracy,  but  RM  model  leads  to  very  poor  results  in  the 
boundary  area. 

Let  us  mention  (see  [35] )  that  theoretical  strength  of  the  boundary  layer 
for  the  model  (1,1,2)  is  -  13.1  which  is  in  very  good  agreement  with 

the  data  in  Table  5.7.  For  additional  details  see  also  [36], 

Summarizing  we  see  that  no  model  guarantees  good  results  and  hence  only 
adaptive  approaches  with  a-posteriori  estimates  could  lead  to  the  accuracy 
with  prescribed  accuracy. 
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The  solution  of  the  three  dimensional  plate  formulation  has  singular 


behavior  in  the  neighborhood  of  the  edges  and  vertices  and  the  exact  character 

of  this  singular  behavior  is  known  [37]  [38].  Although  the  stresses  are  in 

general  singular  in  the  neighborhood  of  all  edges,  the  moments  and  shear 

forces  computed  from  the  three  dimensional  solutions  are  singular  only  in 

the  neighborhood  of  the  corners  of  the  boundary  T  of  u.  This  behavior 

is  also  influenced  by  the  boundary  layer.  Nevertheless  in  the  very  small 

neighborhood  of  the  corner  of  T  (for  simplicity  placed  in  the  origin)  the 

moments  M,  and  shear  forces  Q„  can  be  written  (for  red)  in  the  form 
*  J 

Mjj  =  C^r*1  (0)  +  smoother  terms 

(6.1) 

Q3j  *  Djr*2  +  smoother  terms 

where  r,  0  are  polar  coordinates  and  C^,  Dj  are  analogs  to  the  stress 

intensity  factors  of  the  plane  elasticity. 

Any  model,  for  example  RM,  n-model  or  t  '  Kirchhoff  model  (which  will  be 

addressed  in  Section  7),  yields  the  solution  behavior  of  the  same  form  as  in 

(6.1)  but  the  coefficient  A. ,  i  =  1,2,  and  the  functions  <p.  .(0),  ^.(0) 

*  i » J  J 

depend  on  the  model. 

For  the  analysis  and  the  determination  of  we  refer  to  [26]  where  we 
have  proven  the  following  theorem. 

Theorem  5. 1  i)  For  any  model  (n,n,n+l),  n£l  the  coefficients  in  (6.1) 
are  the  same  as  for  the  3  dimensional  solution. 

ii)  The  Riessner-Mindlin,  as  well  as  the  Kirchhoff  model,  yield  in 
general  the  singular  behavior  which  is  different  when  compared  with  the  3 
dimensional  solution. 
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In  [26]  (see  also  [36]),  we  address  these  singular  behaviors  in  details. 
Here  we  will  Illustrate  the  results  on  only  one  example. 

Consider  the  square  plate  uniformly  bounded,  which  is  clamped  on  the 
bounded  sides  T 2,  T4  and  free  on  the  vertical  sides  1^,  T3  (see  (2.2)). 

The  singularity  coefficients  X^,  X^  for  the  Riessner-Mindlin,  (1,1,2), 
and  Kirchhoff  model  are  given  in  Table  6.1. 

Table  6.1  The  singularity  coefficients  X^,  X^  for  various  models 


RM 

3  DIM 
(1,1,2) 

K 

*1 

0.7583 

0.7112 

1 . 0687 

+i  0.4386 

*2 

1 

1 

0.687 

+i  0.4386 

In  the  case  when  in  (6.1)  is  an  integer,  the  notion  of  smoother 

terms  has  to  be  properly  understood.  Nevertheless,  we  will  not  go  here  in 
details,  although  this  case  is  present  in  our  example.  If  the  coefficient  X 
is  a  complex,  then  (6.1)  has  to  be  understood  as  the  real  part  of  this 
expression  and  the  moments  and  shear  forces  oscilate. 

In  Fig.  6.1,  we  show  in  log-log  scale,  the  moment  on  the  line 

Xj  *  x 2  as  the  function  of  r  (distance  from  the  vertex  of  the  plate).  We 
used  d  *  0.01  and  v  *  0.3.  In  the  figure  we  also  show  the  theoretical 
slope  based  on  the  data  shown  In  the  table  6.1. 
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Figure  6.1  The  moment  M 22  for  the  3  dimensional  solution  and  RM  model. 

In  Figure  6.2  we  show  the  moments  M^,  M ^  as  well  shear  force  Q^2 
for  the  RM  model. 
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Returning  to  the  table  6. 1  we  see  that  the  Klrchhoff  model  yields 
oscilating  behavior  of  the  moment  ^22'  a^so  mentioned  that  for  d — >0, 
all  the  models  converge  to  the  Klrchhoff  solution  Inside  the  domain  w.  Hence 
we  can  expect  that  on  the  line  x^  =  x^  we  will  see  Increased  number  of 
oscillations  in  the  moment  M ^  as  d — >0.  Figure  6.3  shows  the  behavior  of 
the  moment  comPu^ec*  from  3  dimensional  solution  for  d  =  0.0.1.  We  see 
clearly  the  effect  of  the  predicted  oscillation  although  in  a  weak  form. 


Figure  6.3  The  behavior  of  the  moment  for  the  3  dimensional  solution. 

7)  The  Klrchhoff  model. 

The  Klrchhoff  model  describes  the  limiting  solution  of  the  plate 
problem  as  d — >0.  It  can  describe  only  the  restricted  set  of  the  boundary 
conditions.  Let  us  address  this  condition  in  relation  to  the  problem  of  simply 
supported  plate  and  one  would  like  to  know  whether  the  Klrchhoff  model  appro¬ 
ximates  the  hard  or  soft  simple  support.  Because  of  physical  interpretation 
of  the  Klrchhoff  hypotheses,  we  can  expect  that  it  approximates  the  hard 
simple  support.  This  can  be  further  supported  by  the  following  observations. 
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The  Kirchhoff  model  leads  to  the  polygon  paradox  as  the  hard  simple 
support  in  3  dimensional  setting  does. 

Table  7.1  shows  (d  *  0.01,  v  =  0.3)  the  comparison  between  the  Kirchhoff 
and  3  dimensional  formulation 

Table  7.1  The  comparison  of  the  Kirchhoff  and  3  dimensional  formulation. 


Q31 (0.4,0. 5) 

Q31(0.S-0.0.5) 

M12 (0. 5-0, 0. 5) 

Hard  support 

0 

0 

0.0325 

Soft  support 

-9.72 

-6.79 

0 

Kirchhoff 

0 

0 

0.0325 

As  usually  we  denoted 

Q31(0.5-0,  0.5)  =  lim  Q  (x,  0.5) 
x*0.5 

Often  in  the  Kirchhoff  theory  (see  [39]),  the  reaction,  for  example 
V2  =  Q32  (for  x2  =  0.5),  is  computed  from  the  formula 

V  =  fQ  -  aJh£\ 

2  1^32  dxi  ) 

The  formula  tries  to  simulate  the  reaction  for  the  soft  support  by  the 
solution  for  the  hard  support.  Although  this  simulation  is  partly  successful, 
it  does  not  lead  to  the  limiting  value  for  d — >0.  Further  the  approxi¬ 
mation  values  for  Q32(x1>0.5)  cannot  be  predicted  at  all  (see  Table  7.1). 

The  Kirchhoff  model  in  the  case  of  simple  support,  leads  to  the  negative 
concentrated  reaction  in  the  corner.  It  is  interesting  to  compare  the  value 
of  this  reaction  with  the  integral  of  negative  shear  force  (reaction  for  the 
soft  support.  We  get  R  =  0.0290  compared  with  the  Kirchhoff  value  V  =  0.0325, 
l.e.  10%  difference.  For  more  details  we  refer  to  [35]. 
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8.  The  Optimal  Models 

In  the  Section  3,  we  assumed  the  form  (3.1)  as  the  base  for  the  approxi¬ 
mate  solution  of  the  plate  problem.  In  the  Section  4  we  concluded  that,  from 

the  asymptotic  point  of  view,  <p.  ,(?})  should  be  a  polynomials.  We  can  ask  the 

*  *  J 

question  whether  polynomials  are  optimal  when  d  is  fixed.  The  functions 

9?  .(■»»)  will  be  optimal  if  they  give  the  minimal  error  among  all  possible 
i  »  J 

choices  of  <p.  At))  (e.g.  in  the  energy  norm).  In  [31]  it  has  been  shown  that 
1  »  J 

polynomials  are  not  optimal  for  d>0.  The  optimal  <p,  (tj)  are  essentially 

1  *  J 

of  the  form  shat)  where  a  depend  on  the  class  of  function  under  considera¬ 
tion  and  the  thickness  d.  The  proper  choices  are  of  great  importance  for  the 
laminated  plates.  For  more  details  about  the  hierarchical  models  for  laminated 
plates  we  refer  to  [40]. 


9.  Adaptive  modelling  and  a-posterlori  error  estimates 

We  based  our  models  on  the  form  (4.1).  Obviously  the  term 


(HJ 


can  be 


replaced  by  an  arbitrary  polynomial  of  degree  S  j  in  x^.  We  would  like  to 
have  only  a  minimal  number  of  terms  which  will  contribute  to  the  solution  i.e. 
that  |u^  j|  will  decrease  with  increasing  j.  This  can  be  achieved  by 
selecting  j  *  where  Pj  is  the  Legendre  polynomial.  Then  the 

adaptive  approach,  which  selects  different  n  in  the  different  regions  of  the 
plate,  can  be  governed  by  the  requirement  that  |u,  I  >c,  i  =  1,2,3  in  the 
region,  where  it  is  defined.  The  a-posteriori  error  estimation  can  also  be 
extracted  from  the  analysis  of  the  change  when  n^  is  increased.  This  is  in 
the  same  spirit  as  the  error  estimation  in  the  finite  element  method.  In  fact 
we  can  understand  the  modelling  as  special  case  of  the  finite  element  method 
[33].  Nevertheless  we  will  not  go  here  in  more  details.  For  some  results 
about  adaptive  approach  we  refer  also  to  [41]. 
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10.  Numerical  solution  of  the  plate  models. 

Using  the  form  (4.1)  In  the  variational  principle  we  can  understand  it 
as  the  p-version  of  the  finite  element  solution  when  degree  q  Is  used  in 
variable  [33].  In  xj»x2  d°®aln  «  is  partitioned  into  elements  and  the 
h  or  p  or  hp  version  of  finite  element  can  be  used. 

The  discretions  in  the  variables  x^.x^  leads  to  a  locking  problem.  For 
general  theory  of  locking  and  quantitative  assessment  of  the  locking  we  refer 
to  [42].  The  locking  expresses  the  fact  that  as  d — >0,  the  performance 
deteriorates.  Nevertheless  it  is  necessary  to  realize  that,  as  d — >0,  the 
boundary  layer  becomes  stronger  and  hence  it  also  influences  the  performance. 
Hence  in  the  locking  analysis  and  the  numerical  experimentation  the  smoothness 
of  the  exact  solution  has  to  be  independent  of  d. 

Let  us  address  the  solution  of  the  RM  model  (the  (1,1,0)  model  with 

B  =  R).  We  assume  that  in  the  variables  (Xj^)  the  elements  are  rectangular 

of  degree  p  (the  serendipity  elements)  on  square  mesh  of  the  size  h.  Then 

the  rate  of  convergence  due  to  locking  is  hP  *  (measured  in  the  energy  norm), 

while  the  best  possible  approximation  is  of  the  order  h  .  Hence  we  have  a 

loss  of  "one  unit"  in  the  rate.  Nevertheless,  this  "loss"  for  higher  order 

elements  is  visible  only  for  d  very  small.  Let  us  show  an  example.  We 

consider  the  unit  square  plate  n  with  hard  simple  support  (because  the 

s 

boundary  layer  is  very  week  (w  =  0.3)).  Figure  10.1  shows  the  accuracy  of  the 
method  when  uniform  mesh  of  different  sizes  h  is  used.  The  rate  hp  is  also 
shown  in  the  Figure  10.1.  We  see  that  for  d  *  0.01  and  p  =  1  convergence 
is  virtually  not  present  while  for  p  -  4  we  do  not  see  any  loss  of  convergence 
rate  for  all  0.12d£0.01. 
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Hence  the  locking  is  essentially  not  Influenced  by  the  selection  of  the  model. 

It  is  necessarily  to  realize  that  higher  models  show  stronger  boundary 
layer  and  its  relation  to  the  mesh  size  can  lead  to  the  wrong  conclusion 
that  the  locking  disappears  as  d — >0.  Figure  10.3  shows  the  effect  for 
the  model  (3,3,4).  For  more  se  [33]. 


Figure  10.3  The  convergence  of  the  h-version  for  the  (3,3,4)  model  (p=4). 


28 


References 


[1]  Kirchhoff,  G.  [1850]:  Ober  das  Gleichgewicht  und  die  Bewegung  einer 
elastischen  Scheibe  J.  Reine  angew.  Math.  40,  51-58. 

[2]  Ambartsumjan,  S.A.  [1987]:  Theory  of  anisotropic  plates,  Nauka,  Moscow 
(in  Russian). 

[3]  Pane,  V.  [1975]:  Theories  of  elastic  plates.  Noordhoff,  Leyden 
Netherlands. 

[4]  Noor,  A.K. ,  Burton,  W.S.  [1989]:  Assessment  of  shear  deformation 
theories  for  multilayered  composite  plates,  Appl. Mech. Rev. 42,  1-12. 

[5]  Goodier,  J.N.  [1938]:  On  the  problems  of  the  beam  and  plate  in  the 
theory  of  elasticity.  Trans  Roy.  Soc.  Canada  32,  65-88. 

[6]  Friedrichs,  K.O.,  Dressier,  R.F.  [1961]:  A  boundary-layer  theory  for 
elastic  plates.  Comm. Pure  Appl.  Math.  14,  1-33. 

[7]  Goldenveizer,  A.L.  [1962]:  Derivation  of  an  approximate  theory  of 
bending  of  a  plate  by  the  method  of  asymptotic  integration  of  equations 
of  the  theory  of  elasticity.  Prikl.  Math.  Mech.  26,  668-686  (English 
translation  J.  Appl.  Math.  Mech.  (1964),  1000-1025). 

[8]  Reissner,  E.  [1944]:  On  the  theory  of  bending  of  elastic  plates  J. 
Math,  and  Phys.  23,  184-191. 

[9]  Reissner,  E.  [1963]:  On  the  derivation  of  boundary  conditions  for 
plate  theory  Proc.  Roy.  Soc.  Ser.  A276,  178-186. 

[10]  Reissner,  E.  [1985]:  Reflections  on  the  theory  of  elastic  plates.  Appl. 
Math.  Rev.  38,  1453-1464. 

[11]  Naghdl,  P.M.  [1972]:  The  theory  of  shell  and  plates  in  Handbuch  der 
Physics.  Vol.  VIa(2)  (S.  Flugge  and  C.  Truesdell  eds. )  425-640 
Springer  Verlag  Berlin. 

[12]  Clarlet,  P.G.  [1990]:  Plates  and  Junctions  in  elastic  multi-structures. 
An  asymptotic  analysis.  Masson,  Paris,  Springer,  Berlin. 

[13]  Morgenstern,  D.  [1959]:  Herleitung  der  Plattentheorie  aus  der 
dreidinensionalen  Elastlzitatstheorie.  Arch  Rat.  Mech.  Anal. 4, 

145-152. 

[14]  Destuynder,  P.  [1982]:  Comparaison  entre  les  modules  tridimensionnels 
et  bidlmensionnels  de  plagues  en  Alasticite.  RAIR0,  Analyse  Numurique 
15,  331-364. 

[15]  Ciarlet,  P.G. ,  Kesevan,  S.  [1981]:  Two  dimensional  approximation  of 
three  dimensional  eigenvalue  problems  in  plate  theory.  Comp. Meth. Appl. 
Mech.Engrg.  26,  149-172. 


29 


[16]  BabuSka,  I.,  Prager,  M.  [I960]:  Reissner  and  algorithm  in  the  theory 
of  elasticity.  Bull.  Acad.  Polon.  Sci.  Ser.  Sci.  Math.  Astr.  Phys.8, 
411-417. 

[17]  Mindlin,  R.D.  [1951]:  Influence  of  rotatory  inertia  and  shear  on 
flexural  motions  of  isotropic  elastic  plates.  J.  Appl.  Math.  18, 

31-38. 

[18]  Babuika,  I.  [1961]:  The  stability  of  the  domain  of  definition  with 
respect  to  basic  problems  of  the  theory  of  partial  differential 
equations  especially  with  respect  to  the  theory  of  elasticity  I,  II, 
Czechoslovak.  Math.  J.  76-105,  165-203  (in  Russian). 

[19]  BabuSka ,  I.  [1962]:  The  theory  of  small  change  in  the  domain  of 
definition  in  the  theory  of  partial  differential  equations  and  its 
applications.  Proceedings  of  the  Conference  Equadiff.  Prague  1962, 
13-26,  Academia. 

[20]  BabuSka,  I.  [1962]:  The  stability  of  domain  and  the  question  of  the 
formulation  of  plate  problems,  Appl.  Mat.  463-467  (in  German). 

[21]  BabuSka,  I.  Pitkaranta,  J.  [1990]:  The  plate  paradox  for  hard  and  soft 
simple  support  SIAM  J.  Math.  Anal.  21,  551-571. 

[22]  Mazya,  V.G. ,  Nazarov,  S.A.  [1984]:  About  the  Sapondzhyn-Babu§ka  Paradox 
in  the  plate  theory,  Dokl  Akad  Nauk  Arm.  Rep.  78,  127-130  (in  Russian). 

[23]  Arnold,  D.N. ,  Falk,  R.S.  [1989]:  Edge  effects  in  the  Reissner-Mindlin 
Plate  theory.  In  analytical  and  computational  models  of  shells  ed. 

A.K.  Noor  T.  Belitsko,  S.  Simo,  ASME,  71-90. 

[24]  Arnold,  D.N.,  Falk,  R.S.  [1990]:  The  boundary  layer  for  the  Reissner- 
Mindlin  plate  model:  Soft  simply  supported  soft  clamped  and  free  plates 
SIAM  J.  Math.  Anal.  21,  281-312. 

[25]  Haggblad,  B. ,  Bathe,  K.J.  [1991]:  Specifications  of  Boundary  Conditions 
for  Reissner  Mindlin  Theory  Based  on  Plate  Bendings  Finite  Elements. 
Int.  J.  for  Num.Math.  Engrg.  to  appear. 

[26]  Li,  L. ,  BabuSka,  I.  [1991]:  On  the  corner  singularities  of  the 
solutions  of  various  plate  bending  models.  Tech.  Note  BN  1118, 
Institute  for  Physical  Science  and  Technology,  University  of  Maryland, 
College  Park,  Maryland  20742. 

[27]  Burton,  W.S. ,  Sinclair,  G.B.  [1986]:  On  the  singularities  in  Reissner’ s 
theory  for  the  bending  of  elastic  plates.  J.  Appl.  Mech.  53,  220-222. 

[28]  Babuika,  I.  [1990]:  The  problem  of  modelling  the  elastomechanics  in 
engineering.  Comp.  Math.  Appl.  Mech.  Engrg.  82,  155-182. 

[29]  Vogelius,  M. ,  BabuSka,  I.  [1981]:  On  the  dimensional  reduction  method 
I.  The  optimal  selection  of  basis  functions.  Math.  Comp.  37,  31-45. 


30 


[30]  Vogelius,  M. ,  BabuSka,  I.  [1981]:  On  the  dimensional  reduction  method 
II  Some  approximation-theoretic  results.  Math. Comp.  37,  47-68. 

[31]  Schwab,  Ch.  [1939]:  The  dimensional  reduction  method.  Ph.D.  Thesis 
University  of  Maryland,  College  Park,  Maryland  29742,  USA. 

[32]  Szabo,  B.A. ,  Babus>ka,  I.  [1991]:  Finite  element  analysis.  J.  Wiley 
1991. 

[33]  BabuSka,  I.,  Li,  L.  [1990].  The  h-p  version  of  the  Finite  Element 
Method  in  the  Plate  Modelling  Problem,  Tech.  Note  BN-1115,  Institute 
for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742. 

[34]  Cowper,  G.  R.  [1966).  The  shear  coefficient  in  Timoshenko  beam  theory. 
Jour,  of  Appl.Mech.  335-340. 

[35]  Li,  L. ,  Babuska,  I.  [1990].  On  the  Riessner-Mindlin  and  the  (1,1,2) 
model  for  the  clamped-in  plate  problem.  Tech.  Note  BN  1117,  Institute 
for  Physical  Science  and  Technology.  University  of  Maryland,  College 
Park,  Maryland  20742. 

[36]  BabuSka,  I.,  Li,  L.  [1990].  The  problem  of  plate  modelling. 

Theoretical  and  computational  results.  Tech.  Note  BN  1116,  Institute 
for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742. 

[37]  Dauge,  M.  [1988]:  Elliptic  boundary  value  problems  on  corner  domains. 
Lecture  Notes  in  Math.  1341,  Springer,  New  York  1988. 

[38]  Ptersdorff,  T.  [1980]:  Randwertprobleme  der  Elastizitatstheorie  fur 
Polyder  Singularitaten  und  Approximation  mit  Randenelement  Methoden, 
Ph.D.  thesis,  T.H.  Darmsdadt,  FRG. 

[39]  Timoshenko,  S.  Woinowski-Krieger,  S.  [1959]:  Theory  of  plates  and 
shells,  Sec.  ed.  McGrow-Hill. 

[40]  Babuika,  I.,  Szabo,  B.A.,  Actls,  R.L.  [1990]:  Hierarchic  Models  for 
laminated  composites  Rep.  WO/CCM-90/40,  Center  for  Computational 
Mechanics,  Washington  University,  St.  Louis  M0. 

[41]  Vogelius,  M. ,  Babu§ka,  I.  [1981]:  On  a  dimensional  reduction  method 
III.  A-posterior  error  estimation  and  an  adaptive  approach.  Math. 

Comp.  37,  156,  361-383. 

[42]  BabuSka,  I,  Suri,  M.  [1990]:  On  locking  and  robustness  in  the  finite 
element  method.  Tech.  Note  BN  1112,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland. 


31 


The  Laboratory  for  Numerical  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Maryland,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals: 

•  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 

in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology.  '  * 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  Babuska,  Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742. 


